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We present a study of the electronic ultrasound attenuation a in the unconventional superconductor 
Sr2Ru04. The power law behavior of a at temperatures down to T c /30 clearly indicates the presence 
of nodes in the gap. In the normal state, we find an enormous anisotropy of a in the basal plane 
of the tetragonal structure. In the superconducting state, the temperature dependence of a also 
exhibits significant anisotropy. We discuss these results in relation to possible gap functions. 
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The symmetry of the order parameter A is of funda- 
mental importance in understanding the nature of su- 
perconductivity. To this day, only two symmetries have 
been unambiguously identified amongst all known super- 
conductors: s-wave symmetry in conventional supercon- 
ductors and d-wave symmetry in some of the high-T c 
cuprates. Both cases involve spin-singlet Cooper pairing. 
Heavy-fermion and organic superconductors are widely 
believed to be unconventional (i.e. not s-wave); how- 
ever, for a variety of technical reasons, the symmetry 
of A has not been firmly established in any of these sys- 
tems. In this context, Sr2Ru04 has emerged as a highly 
promising candidate for the study of unconventional su- 
perconductivity jjj because of its thoroughly character- 
ized Fermi Surface (FS) 0] and the availability of sizable 
high-quality crystals. Moreover, it may well prove to be 
the first clear example of a spin-triplet superconductor 

I- 

The most direct evidence for spin-triplet pairing is the 
absence of detectable change in the Pauli susceptibil- 
ity below T c as measured by 17 Knight shift jl and 
neutron scattering H. In addition, the appearance of a 
small random spontaneous internal field in the supercon- 
ducting state H indicates that A breaks time-reversal 
symmetry, a conclusion supported by analysis of the 
square flux lattice 0. The simplest A consistent with 
these two properties possesses an orbital wave function 
d = z(k x + iky). This p-wave state has an isotropic gap 
(without nodes), and for several years it was viewed as 
the leading candidate for the superconducting order pa- 
rameter of Sr2Ru04 S. 

However, in the last year it has become clear that an 
isotropic gap is inconsistent with the power law T depen- 
dence of several quantities observed deep in the super- 
conducting state |)-12|. Such behavior is most naturally 
explained by line nodes in the gap. Thus the detailed 



symmetry of the order parameter in Sr2Ru04 is still very 
much an open question. In order to determine the orbital 
symmetry of A, one requires an angle-resolved probe to 
reveal the location of nodes in the gap on the FS. One 
approach is to seek an anisotropic response to an ap- 
plied magnetic field. Such an effect has been observed in 
heat transport with the field rotating in the basal plane. 
However, the anisotropy was quite small and was deemed 
incompatible with vertical line nodes pj. In the study 
presented in this Letter we use the intrinsic angular res- 
olution of ultrasound attenuation to directly probe the 
anisotropy of the quasiparticle excitation spectrum of su- 
perconducting Sr2Ru04 in zero magnetic field. 

Moreno and Coleman have emphasized the power of 
ultrasound attenuation in studies of anisotropic super- 
conductors |0|]. In their model calculations, performed 
in the hydrodynamic limit where the electron mean free 
path I is much shorter than the sound wavelength A, the 
low T power law behavior of a{T) is shown to depend 
strongly on the direction of sound propagation q and po- 
larization e relative to the nodes. For example, in the 
case of a 2D gap with k\ — ky symmetry and for trans- 
verse sound with q and e in the basal plane, quasipar- 
ticles at nodes along the (110) are "inactive" when q || 
(110), yielding a strong T 3,5 power law, whereas they are 
maximally "active" in attenuating sound for q || (100), 
resulting in a weak T dependence. In this way ultra- 
sound attenuation experiments can locate nodes in the 
gap, as was done for the heavy-fermion superconductor 
UPt 3 @. 

We have performed measurements of longitudinal and 
transverse ultrasound attenuation in Sr2Ru04 in the 
temperature range 0.04 - 4 K, with sound propagating 
along the (100) and (110) directions in the basal plane of 
the tetragonal crystal structure. Our experiments were 
carried out on oriented pieces cut from a single high- 
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quality crystal grown by the travelling solvent float zone 
technique . T c defined as the peak of the magnetic 
susceptibility (x") is 1.37 K. Samples were polished to 1 
fim roughness, with two opposite faces whose parallelism 
was estimated to be better than 1.5 /mi/mm. The align- 
ment of the polished faces relative to the crystal axes 
was determined by Laue back reflection to be less than 
0.5° off axis. (For propagation along (100), we estimate 
the misalignment to be 0.2-0.3°.) The crystals were 3.98 
and 4.23 mm long, and had overlapping regions of the 
two parallel faces of ~ 2 mm 2 . Pulse-echo ultrasound 
was performed using a homebuilt spectrometer and com- 
mercial (30 MHz) LiNbC>3 transducers, bonded to the 
crystals with a thin layer of grease. With this technique, 
the absolute sound velocities can be measured to a few 
percent, but only a relative measurement of the electronic 
attenuation is possible because one cannot distinguish 
power losses from other origins. We denote the modes as 
follows: T100, T110, L100, and L110, where the polar- 
ization (T for transverse, L for longitudinal) is always in 
the plane. The indices denote the propagation direction. 
From the echo spacing, the sound velocities were found to 
be 3.30, 2.94, 6.28, and 6.41 mm/^s, respectively, within 
3%. (These are somewhat different from previous reports 
p6| .) Typically at low temperature for both the T100 
and LI 10 modes we recorded more than 50 echoes, while 
L100 had at least 20 and T110 at least 2. Preliminary 
measurements were accomplished in a 3 He cryostat and 
were consistent with data taken in a dilution refrigerator 
up to 4 K. The results were also independent of excitation 
power. 

An important parameter that controls the behavior of 
the attenuation is the ratio X/l. Theoretical approaches 
are quite different depending on the value of this param- 
eter. The electronic mean free path below T c is on the 
order of 1 /zm, as estimated from either T c or resistiv- 
ity fl7|] . From this we expect that the crossover region 
ql ~ 1 (where q = 2tt/X is the sound wavevector) from 
the hydrodynamic to the quantum regime occurs at 0.5 
(1.0) GHz for the transverse (longitudinal) modes. The 
theoretical dependence of a on frequency is v 2 in the hy- 
drodynamic regime and v in the quantum regime. Fig. [j] 
shows data for two modes at two frequencies, normalized 
at T c . The values of a(T c ) obey approximate v 2 scaling, 
and the shape of a{T) is independent of v and repro- 
ducible at many frequencies up to 500 MHz. We conclude 
that all our data are in the hydrodynamic limit. 

In this limit, the attenuation is related to the 4 th rank 
electronic viscosity tensor (77) via the relation |Tj| : 

(2^) 2 _ . (2nv) 2 _ 

a = ~—Vijkieiqje k qi = ~— ?7modo, (1) 

PC A S pc A s 

where p is the mass density and c s is the sound velocity 
for the mode with polarization and propagation direc- 
tions e and q. This defines ^ mo dc as the linear combi- 
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FIG. 1. Ultrasound attenuation a(T) normalized at T c , 
both in the superconducting state (B = 0) and in the nor- 
mal state (B = 1.5 T 1 c). We show both T100 and L100 
at different frequencies to demonstrate that the shape of the 
curves does not depend on frequency. 

nation of rjiju components for a particular mode. Un- 
der tetragonal symmetry there are 6 independent com- 
ponents of 77, but only three determine the attenuation of 
the four in-plane modes, so that rymodc is given by 771212, 
771111, § (171111 ~ 771122), and §(771111 + 771122 + 2771212) for 
T100, L100, T110, and L110, respectively. We present 
TTmodc in Fig. ^. The curves show that L100 and T110 
have the same shape and almost the same amplitude. 
The LI 10 mode has a similar shape but a smaller am- 
plitude. However, the T dependence of the T100 mode 
is clearly weaker and the magnitude of 77x100 is about 
1000 times smaller than 77L100 or 77x110- We note that 
the alignment of the crystal faces mentioned above is 
sufficient to yield the intrinsic 77x100 j[9)- Such a dra- 
matic anisotropy of the normal state electronic 77 is un- 
precedented. Next we consider possible origins for this 
anisotropy. 

The viscosity is related to the deformation of the FS 
caused by the sound wave's strain field, weighted by the 
electronic relaxation time r pC( |. The anisotropy of 77 
may in part come from angular dependence of r or the 
topology of the FS which is composed of three separate 
sheets (labelled a, (3, and 7) ||. Although the 7 sheet 
is fairly isotropic, the a and (5 sheets have flat regions 
that could make significantly different contributions to 
the different components of ?;. Intcrband effects at the 



2 




0.0 0.5 1.0 1.5 2.0 



Temperature (K) 

FIG. 2. Viscosity for the four in-plane modes. The viscos- 
ity r\ is related to the attenuation by r\ — apct. / '{2nv) 2 , using 
a density p — 5.9 g/cm 3 and the sound velocities (c s ) given 
in the text (poise = g/cm s). The data is extrapolated to K 
to obtain the zero of the vertical axis. 

point of near contact between the three sheets may also 
be important. However, the anisotropy is so large that 
it probably requires significant variation in the electron- 
acoustic phonon coupling for the different modes, possi- 
bly originating in the highly directional nature of the Ru 
Ad I O 2p hybrid orbitals that compose the conduction 
bands. 

To facilitate discussion of a(T) below T c , we continue 
with some remarks about the observed viscosities. The- 
ory of 2D electronic viscosity j2l]] suggests a close re- 
lation between two of the relevant components, namely 
771m = — 771122 for certain simple FS. The data pre- 
sented in the top panel of Fig. || shows that this rela- 
tion is nearly fulfilled in S^RuCU, in both magnitude 
and T dependence. The slight imperfection of this sym- 
metry should not be disturbing, since, e.g., the real FS 



is not exactly 2D. This relation also implies that 7y mo de 
in the bottom two panels should be nearly equal. The 
apparent dramatic violation of this prediction is, how- 
ever, a trivial consequence of the slight imperfection of 
the symmetry together with the enormous anisotropy 
(1^71212 1 < himl, I»7ii22 1 ) - While r7 L n contains all three 
viscosities, it is dominated by the slight inequality be- 
tween 77L100 and 77T110 (as the difference is still much 
larger than 771212). Because the observed fjLno(T) is so 
similar to 77lioo(T), this difference is primarily one of 
magnitude. Thus there are effectively only two T depen- 
dences for the in-plane viscosities, one represented exclu- 
sively by 77x100 and the other most clearly by ttlioo- We 
compare these in detail next. 
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FIG. 3. ct(T) divided by the normal state attenuation 
vs reduced temperature. The solid line is the simple BCS 
prediction. Inset: data up to T c , compared with calculations 
for /-wave discussed in the text [El). 

The attenuation for these two modes is shown in Fig. ||, 
normalized by their normal state curves a^r{T) measured 
in a magnetic field (1.5 T, approximately in-plane) ex- 
ceeding the upper critical field. A simple BCS model 
(weak-coupling isotropic s-wave gap) or a clean fully 
gapped p-wave state |22j give an exponential decrease 
of a(T) at low T, such that below ~ T c /4 the atten- 
uation is essentially zero (solid line in Fig. ||). This is 
obviously not the case here, where we see power laws 
persisting down to T c /30, providing compelling evidence 
for the presence of nodes in the gap. It is also clear that 
there is anisotropy in the T dependence of the viscosity 
in the basal plane. The weakest T dependence occurs 
for the T100 mode. Below 0.7 K, the normalized ?7tioo 
and t7lioo are well described by the power laws 0.53i 14 
and 0.38t 18 respectively, where t = T/T c . The latter 
power law is consistent with another recent measurement 
of t7tho p3| ]. The anisotropy revealed in Fig. || is much 
larger than any other kind of anisotropy observed so far 
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in the basal plane of superconducting Sr2Ru04. Next we 
discuss the anisotropy of a(T) in relation to the nodal 
structure of the gap. 

At the outset we note that the huge anisotropy in the 
normal state viscosity is not included in theories forming 
the basis of existing calculations of a(T) in the super- 
conducting state. Thus, we can only draw tentative con- 
clusions from a comparison of these theories to the data. 
We identify two types of order parameter that have been 
proposed: those with vertical line nodes (e.g. pH) and 
those with horizontal line nodes (e.g. p4|,E5|]). We now 
consider these in turn. 

If the angular sensitivity articulated by Moreno and 
Coleman survives in a theory correctly accounting for 
the normal state viscosity then the anisotropy of Fig. || 
is naturally interpreted in terms of vertical nodal lines. 
In this case, the observed anisotropy of the T dependence 
between L100 and T100 is qualitatively consistent with a 
A of k% — ky orbital symmetry in 2D such as found in the 
cuprates. However, given that a d-w&ve order parameter 
is inconsistent with the combination of spin-triplet sym- 
metry and time-reversal symmetry breaking observed in 
Sr2Ru04, one needs to invoke accidental nodes, such as 
those of an /-wave state with d = z(k x + ik y )(k\ — ky) 
p4j . In the inset of Fig. ||, we compare our normalized 
data for L100 and T100 with the calculation of Graf and 
Balatsk y for such an /-wave gap on a single isotropic 
2D FS |2(J. In detail the T dependences differ signifi- 
cantly from the simple model predictions. There is good 
overall agreement for the L100 mode; however, the low 
T power law exponent is much lower than expected (1.8 
vs 3.5). The data for the T100 mode is always below the 
calculated curve, but, in contrast, the exponent is close 
to the predicted 1.5 (T^j2l[]. These significant differences 
may be the consequence of the real FS of Sr2Ru04 or of a 
more complex nodal structure. Detailed calculations are 
needed to evaluate how strongly the observed anisotropic 
T dependence favors such a A. 

For a simple horizontal line node, one would not ex- 
pect anisotropy in the plane. Thus we must presume a 
separate source of the observed anisotropy. One possibil- 
ity is a (non-nodal) 4-fold variation of |A| in the plane. 
The magnitude of the modulation required to explain 
our data could be computed and the resulting gap func- 
tion compared with other measurements such as the elec- 
tronic specific heat. However, at sufficiently low T, we 
still expect the power law exponent to be isotropic in this 
scenario. 

In summary, the power law T dependence of the ultra- 
sound attenuation persisting down to T c /30 establishes 
unambiguously the existence of nodes in the gap func- 
tion of Sr 2 Ru04. These power laws and the anisotropy 
of the attenuation in the superconducting state excludes 
the possibility of a pure p-wave state, but may prove to be 
consistent with a state that has either vertical nodes on 
the diagonal or a horizontal line of nodes in conjunction 



with significant 4-fold gap modulation. A detailed cal- 
culation taking into account the real (multi-sheet) FS of 
Sr2Ru04 and the large normal state anisotropy is needed 
to properly extract the full order-parameter information 
contained in the anisotropy and temperature dependence 
of a(T). 
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Note in Ref. |2l| the /-wave state is d = z(k x +ik y )(k x k y 



which is equivalent to z(k x + ik y )(k x - 
about z. 
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